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Density-functional calculations are used to verify the atomic structure of the hexagonal In/Si(lll)- 
{V7 X %/3) surface, which has been considered to represent an ultimate two-dimensional (2D) limit 
of metallic In overlayers. Contrary to the prevailing assumption, this surface consists of not a single 
layer but a double layer of In atoms, which corresponds to a hexagonal deformation of the well- 
established rectangular In double layer formed on Si(lll)-(%/7 x y/3) [Park et al., Phys. Rev. Lett. 
109, 166102 (2012)]. The same double-layer thickness accounts well for the typical coexistence of the 
hexagonal and rectangular phases and their similar 2D electronic structures. It is thus conclusive 
that, regardless of rectangular or hexagonal, the In/Si(lll)-(%/7 x %/3) surface does not represent 
an one-atom-thick In overlayer. 

PACS numbers: 68.43.Fg, 68.47.Fg, 73.20.At 


One-atom-thick metal overlayers grown on semicon¬ 
ductor surfaces have attracted great attention as model 
systems ideal for e^loring intriguing low-dimensional 
metallic properties One representative system of 

recent interest is the In/Si(lll)-(\/7x \/3) surface, where 
the In overlayer was generally assumed to be one atom 
thick and so represent an ideal two-dimensional 

(2D) limit of metallic In properties. Fascinating elec¬ 
tronic features of the In overlayer, including a nearly- 
free-electron Fermi surface, Q, superconducting transi¬ 
tions ^ 111, and an intriguing metallic transport behav¬ 


ior 


12|, have been explored and referred to as revealing 


the ultimate 2D limit. With little structural information 
about the In over layer, however, its actual layer thickness 
has long been an open question. 

Particularly interesting in this regard is a double-layer 
picture for the In/Si(lll)-(-\/7 x \/3) surface, contrary 
to the prevailing single-layer assumption. In a recent 
density-functional theory (DFT) study [l^, of the two 
different In-derived (\/7 x v^) phases distinguishable in 
scanning tunneling microscopy (STM) experiments [i- 
[^, the more representative, rectangular phase (hereafter, 
•\/7-rect) was verified to contain an In double layer [see 
Fig. 1(a)]. Consequently, the remaining hexagonal phase 
(•\/7-hex) has been spotlighted as a single-layer alterna¬ 
tive, because it typically developes a little earlier than 
the v^-rect phase when prepared by In draosition on 
the Si(III)-(7x7) surface at about 400 °C [3,0. Since 
the \/7-hex phase evolves between the 4x1 and \/7-rect 
phases in the coverage-dependent growth process UM, 
its In cov erag e would be in between 1.0 ML of the 4x1 
phase 15, 16] and 2.4 ML of the v^-rect phase [10 (here, 
1 ML refers to one In atom per surface Si), but the In 
coverage and thickness is yet to resolve. 

In early STM studies, Kraft et al. suggested a 1.0 ML 
single-layer model for the \/7-hex phase [tIIsI, where the 
observed five STM protrusions per (\/7 x v^) unit cell 
were attributed to five In atoms (corresponding to 1 ML). 
This experimental model was more quantified by Shang 


et al. in a recent DFT study [17|. Another single-layer 
model with an In coverage of 1.2 ML was lately proposed 
by Uchida et al. [l0, based on a similarity of their DFT 
simulations to the experimental STM image [ 0 . The 
In coverage itself is also a matter of debate as seen in 
two recent experiments: while Uchihashi et al. [Tol ] sup¬ 
ported the 1.2 ML single-layer model of Uchida et al. 
|18| based on the similarity of STM images, Yamada et 
al. 10 estimated as about 2 ML in their reflection high- 
energy electron diffraction (RHEED) study. Any single¬ 
layer model, however, has an inherent experimental prob¬ 
lem: in their STM study, Kraft et al. observed that the 
coexisting \/7-hex and V^-rect phases have almost the 
same STM heights (with a small difference by 0.1-0.2 A) 
and also show very similar electron tunneling spectra [0 . 
Previously, these similarities were argued to ensure that 
both phases have a single In layer [0 , but now the recent 
double-layer verification of the v^-rect phase [10 should 
imply a double-layer model for the v^-hex phase as well. 

Here, we propose a new structural model for the hexag¬ 
onal In/Si(lll)-(\/7 X \/3) surface. This model features 
a double layer of 2.4 ML In atoms, which actually corre¬ 
sponds to a hexagonal deformation of the In double layer 
of the v^-rect phase. In what follows, the soundness 
of the hexagonal double-layer model is demonstrated by 
DFT calculations. 

We perform DFT calculations using the Vienna ab- 
initio simulation package within the generalized gradient 
approximation [20j and the projector augmented wave 


method 21 


|. The Si(lll) surface is modeled by a 
periodic slab geometry with six atomic layers and a vac¬ 
uum spacing of about 12 A. The calculated value 2.370 
A is used as the bulk Si-Si bond length. Indium atoms 
are adsorbed on the top of the slab, and the bottom is 
passivated by H atoms. The electronic wave functions 
are expanded in a plane-wave basis with a kinetic energy 
cutoff of 246 eV. A 4x6x1 fc-point mesh is used for the 
X v^) Brillouin-zone integrations. All atoms but the 
bottom two Si layers are relaxed until the residual force 
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FIG. 1: (Color online) Double-layer models for In/Si(lll)-('\/7 x v/S)- (a-c) One rectangular (Ref. \M) and two hexagonal 
models. Large (small) balls denote In (Si) atoms. Solid and dashed lines denote (%/7 x %/3) and effective In-(lxl) unit cells, 
respectively. Numbers denote optimized bond lengths and layer spacings. (d) Phase boundary between V^-rect and HEX2. 


components are within O.OI eV/A. 

In light that the i^-hex and v^-rect phases have sim¬ 
ilar STM heights [7|, we examine first the possibility 
of hexagonal deformations of the double-layer \/7-rect 
phase. Figure 1 shows two easy-to-derive hexagonal mod¬ 
els. One is obtained from the \/7-Tect phase by a mono¬ 
clinic deformation of the rectangular In array along the 
y/Z direction, as captured in Fig. 1(b) (hereafter, HEXl). 
The other (HEX2) is obtained from HEXl by a mono¬ 
clinic deformation along the y/l direction as shown in 
Fig. 1(c). In our calculations, while HEX2 is locally sta¬ 
ble, HEXl is unstable and relaxes back to the v^-rect 
phase with no energy barrier. We also examined variants 
of HEXl and HEX2 by allowing lateral displacements of 
In atoms along the or y/l direction, but all of them 
were found to converge to either the v^-rect phase or the 
HEX2 model. Figure 1(c) shows the optimized HEX2 
structure, which forms almost a regular In lattice (ai = 
3.23 A and 02 = 3.41 A): the In-In and In-Si interlayer 
spacings are 2.47 A and 2.61 A in average, respectively, 
comparing well with the values, 2.42 A and 2.60 A, of the 
y/l-rect model. In energetics, the HEX2 model is as sta¬ 
ble as the v^-rect model with a slightly higher formation 
energy by 0.05 eV per [y/l x v^) cell, being a promising 
model for the v^-hex phase. 

Figure 2 shows the simulated STM image of the HEX2 
model, which indeed compares well with the experimen¬ 
tal feature of two bright and three weak spots per unit 
cell In our simulation, two of the three weak spots 
represent the In-In bonding states rather than individual 
In atoms, thus accounting for the five (not six) protru¬ 
sions from six surface In atoms per unit cell. In Fig. 2(c), 
the calculated charge corrugation is in quantitative agree¬ 
ment with the STM line profile : the maximum height 



(a) Experiment (b) Simulation 


(c) Line profile 3.4 A - experiment 



FIG. 2: (Color online) STM cornparison for the HEX2 model, 
(a) STM image taken from Ref. Q]. (b) Simulated image rep¬ 
resenting the surface of constant density with p = 5x10”“* 
e/A^. For better comparison with the low-bias STM image, 
the calculated Fermi level was shifted up by 0.4 eV. The inset 
shows a contour map drawn with a uniform increment of 0.05 
A, where filled circles denote the topmost In atoms, (c) Ex¬ 
perimental (Ref. M) and simulated line profiles taken along 
the arrow marked in (b). 


difference between the peaks (0.23 A) and the interval of 
the two highest peaks (3.32 A) compare well with the 
experimental values of 0.25 A and 3.4 A, respectively. 

It is interesting to compare the topographic heights of 
the HEX2 model and the y/l-rect phase. When simulated 
at a same bias voltage and charge density, the HEX2 
model appears lower by 0.24 A than the y/l-iect model. 
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(b) PDOS 
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FIG. 3: (Color online) (a) DOS for the present \/7-hex 

model, the %/7-rect model of Ref. [T^, and the 1.2 ML model 
of Ref. [T^ . (b) Projected DOS of the In p and s orbitals for 
y/7-hex and V^-rect. 


in good agreement with the STM observation that the 
•\/7-hex phase appears a little lower by 0.1-0.2 A Q. On 
the other hand, the 1.2 ML single-layer model of Uchida 
et al. [l^ was found to appear far lower by 2.21 A. 

Figure 1(d) displays another fascinating feature of 
the HEX2 model: it could form an atomically sharp 
phase boundary with the coexisting \/7-rect phase. This 
boundary matching is achieved by a 60 ° clockwise ro¬ 
tation of the HEX2 model in Fig. 1(c). The resulting 
angular relation matches well the STM observation that, 
when the V^-hex and v^-rect phases coexist, the angle 
difference of their unit cells is 60° Hi. Their typical 
coexistence Q is also explained by the aforementioned 
energetics: since the HEX2 and '/T-rect phases are ener¬ 
getically comparable, both phases would locally develope 
during the In deposition at high temperatures. 

The HEX2 model, consisting of a quasihexagonal In 
double layer, is not only energetically stable but also re¬ 
veals sufficiently sound STM features, so we readily pro¬ 
pose it as the long-sought structural model for the a/T- 
hex phase and further explore its electronic structure. 

Figure 3 shows the calculated density of states (DOS) 
of our v^-hex model. Noticeable is that the '/7-hex 
model reveals almost the same DOS spectrum as the 
'/f-iect model, well reflecting their common structural 
nature (i.e., the same In coverage of 2.4 ML and similar 
double-layer structures), while the 1.2 ML single-layer 
model of Uchida et al. [l^ does not. Our v^-bex model 
is thus compatible with the experimental report that the 
'/f-hex and v^-rect phases reveal similar scanning tun¬ 
neling spectra 0 . There is of course a little difference in 
the DOS curves, reflecting the underlying structural dif¬ 
ference: the small peak at +0.28 eV in the v^-rect model 



FIG. 4: (Golor online) (a) Band structure of %/7-hex. Gircles 
denote the In states that contain more than 50% of charge in 
the In double layer: the states containing more than 40% of 
charge in the top In layer are emphasized by filled circles. 
Solid (green) curves denote the band structure of the (1x1) 
freestanding In double layer. A dashed (red) curve denotes 
a parabolic fitting to the bottom part of the In-derived free- 
electron-like band. Here, the Fermi level was set to zero, (b) 
Fermi contours. The [y/f x v^) and In-(lxl) Brillouin zones 
are denoted by solid and dashed lines, respectively. The arrow 
denotes the line for the band-structure calculation in (a). The 
large (red) circle represents the 2D Fermi circle constructed 
from the parabolic bands described in (a). In the lower panel, 
the (%/? X v^) zone folding of the Fermi circle accounts for 
most of the complex Fermi contours. 


is a little shifted up to +0.53 eV in the v^-hex model. As 
seen in Fig. 3(b), the difference stems from the Px and 
Py orbitals since their lateral interactions are the most 
affected by the in-plane rectangular-to-hexagonal defor¬ 
mation. 

Figure 4 shows the band structure of the v^-hex 
model. It is metallic with a broad distribution of In- 
derived states. Interestingly, the In band structure could 
be traced from those of the (1x1) freestanding In dou¬ 
ble layer [mentioned in Fig. 1(c)]: the (\/7 x v^) zone 
folding of this (1x1) result accounts well for the com¬ 
plicated In bands and the band gap between —3.36 eV 
and —1.75 eV of the y/f-hex phase. A parabolic fitting 
(with an effective mass of 0.92 m/ connects well the In 
states from the bottom at —8.13 eV to over the Fermi 
level. These free-electron-like states were found evenly 
distributed on both of the top and second In layers (not 
shown here), implying that the double-layer unit may be 
essential for the 2D free-electron nature. This surface has 
similar parabolic bands along other k directions, thereby 
leading to a circular Fermi surface: Its radii are 1.39 A 
and 1.36 A along the and ky directions, respectively. 
Figure 4(b) shows the large Fermi circle in the extended 
{V7 x \/3) Brillouin zones, and the actual Fermi contours 
are mostly accounted for by the the {y/f x v^) zone fold¬ 
ing of the large Fermi circle, as seen in the bottom panel. 
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This 2D free-electron nature of the v^-hex phase is sim¬ 
ilar to that of the \/7-rect phase: in our calculations, 
the \/7-rect phase has the band bottom of —8.18 eV, the 
effective mass of 0.95 m^, and the Fermi circle of 1.4 


It is worth mentioning that there is one more (\/7 x 
■\/3) hexagonal phase, differing from the present \/7- 
hex phase mepared at high temperatures (HT) of about 
400 °C [3, Q. Saranin et al. Q obtained a hexagonal 
In/Si(lll)-(-\/7 X v^) phase by room-temperature (RT) 
In deposition on In/Si(lll)-(-\/3 x -x/S) and identified it 
with the HT v^-hex phase on the basis of similar STM 
images. This identification, however, is open to ques¬ 
tions. First, the STM comparison of both phases was not 
fairly done: the RT phase was examined at high biases of 
-1-0.5 V and +2 V j24l|. while the HT phase at as low as - 
0.2 V and -0.12 V [Sj. Indeed, when taken at about +2 V, 
our STM simulations for the v^-hex model do not repro¬ 
duce the STM image of the RT v^-hex surface. Second, 
the RT and HT v^-hex phases appear differently in topo¬ 
graphic height. The HT phase appears higher by about 
2.0 A in STM topograph than the reference (4 x 1) phase 
0 whereas the RT phase shows only a small height differ¬ 
ence of 0.5 A as found in a recent atomic-force microscopy 
image [2^ , which is a strong implication of different \/7- 
hex phases. Finally, the RT phase has a metastable na¬ 
ture: it transforms into a (\/7 x '/j) phase during cooling 
in the range from 265 to 225 K [2^, while the HT phase 
is stable at cryogenic temperatures as observed in super¬ 
conductivity experiments 11 ,lii. Thus, as a metastable 
intermediate phase, the RT phase should be distinguished 
from the stable HT phase. 


In summary, we proposed a double-layer structural 
model for the hexagonal In/Si(lll)-(\/7 x \/3) surface, 
based on the quantitative microscopic and spectroscopic 
examinations by DFT calculations. This double-layer 
thickness with the In coverage of 2.4 ML is the com¬ 
mon structural feature underlying the typical coexistence 
of the hexagonal and rectangular phases and their sim¬ 
ilar electronic structures close to a 2D free-electron gas. 
Therefore, the In/Si(lll)-(\/7 x v^) surface, either rect¬ 
angular or hexagonal, does not represent a single-layer 
limit of metallic In overlayers: such an ultimate 2D limit 
is yet to achieve. 
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